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Abstract 

Motivated by the recent measurement of the like-sign dimuon charge asymmetry, 
we investigate anomalous CP violation in the Bg — Bg mixing within the supersym- 
metry. We show that when gluino diagrams dominate supersymmetry contributions, 
it is very difficult to realize a large Bg — Bg mixing phase under the constraint from 
electric dipole moments barring cancellations. This constraint can be ameliorated 
by supposing superparticles decoupled. In this limit, we find that it is possible to 
achieve the large CP asymmetry, and the branching ratio of Bg l-i-^ fJ'~ tends to 
become sizable. 



1 Introduction 



A lot of experimental efforts have confirmed the violation of the CP symmetry in the 
quark sector such as in the K and B mesons. Although the CP violation is sensitive to 
physics beyond the standard model (SM), most of the CP- violating phenomena are well 
explained within the framework of the SM. In spite of the experimental constraints, a 
room is still left for the h — s mixing sector [1] . 

The Tevatron experiments study possible effects of the CP violation in the B system. 
The DO collaborations recently reported a measurement of the like-sign dimuon charge 
asymmetry. Interpreting the result by the mixing of the neutral B mesons, we obtain 
(A^Jexp = -(9.57±2.51±1.46) x lO'^ [2l[3]. Based on the Tevatron data, the asymmetry 
is also expressed in terms of the semileptonic asymmetry as ~ O.Sa^^ + 0.5a^; [21 [3]. 
Noting that the semileptonic asymmetry is related to the B meson mixing, the dimuon 
charge asymmetry is predicted by the SM as (A^Jsm = (— 2.3lo;g) ^ 10""^ @]- Interestingly, 
the experimental result is larger than the SM prediction at the 3.2(T level. Combined with 
the other experimental data which are sensitive to the B — B mixing, the dimuon charge 
asymmetry favors an anomolous CP phase arising in the B^ meson mixing [5]. 

As a new source of the CP violation, the supersymmetric (SUSY) extension of the SM 
involves soft SUSY breaking parameters. In particular, CP violating phases generally give 
arise in flavor- changing components of the soft scalar masses. They contribute to flavor- 
changing neutral currents (FCNCs) with CP violations. Actually, it has been argued that 
CP- violating h — s squark mixings are a suitable candidate for the anomalous CP violation 
which is favored by the like-sign dimuon charge asymmetry [6l [TJ [8] • 

However, flavor and CP violations have been restricted experimentally. Especially the 
electric dipole moments (EDMs) are sensitive to anomalous CP violations. It was pointed 
out that CP-violating h — s squark mixings induce the chromo electric dipole moment 
(CEDM) of the strange quark and are tightly constrained by the atomic and neutron 
EDMs [9]. Nevertheless, this constraint has been often discarded in previous works on 
CP-violating FCNCs, because estimation of EDMs includes potentially large hadronic 
uncertainties and may also be suppressed by (accidental) cancellations. In contrast, we 
investigate the Bg — Bs mixing phase with seriously taking account of the EDM constraint. 
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We will show that it is very difficult to realize such a large CP asymmetry unless we assume 
cancellations among the SUSY contributions. 

A large CP violation may be reconciled with the EDM bound by supposing cancella- 
tions among SUSY contributions to the EDMs. This may be indeed a possible solution, 
because the phase structure of the SUSY contributions to the Bs — Bs mixing differs from 
that of the EDMs [10]. However, since the EDM bound is too tight, the cancellation is 
quite easily upset by radiative corrections unless they are protected by symmetries. 

Apart from the cancellations, a simple solution to avoid the EDM bound is to decouple 
the (colored) superparticles. Heavy superparticles suppresses most of the SUSY contribu- 
tions to FCNCs as well as those to the EDMs. In this limit, the FCNCs can be induced 
by mediating the heavy Higgs bosons at the BLO level [H]. Actually, non-holomorphic 
(flavor-changing) interactions can be left sizable in the decoupling limit [121 [13]. We will 
show that there is a wide parameter region which realizes the large CP violation in the 
Bs — Bs mixing with alleviating the EDM constraints. In addition, we will predict a 
relatively large branching ratio of the Bs /i/i decay, and it is expected to be measured 
in the LHCb experiment. 

2 CP Phase in Bs Meson Mixing 

The like-sign dimuon charge asymmetry implies an anomalous CP-violating contribution 
to the B meson mixing [21 |3]. The asymmetry is expressed by the i?^ — Bq mixing 
parameters through the semileptonic asymmetry, a^i = Im(r^2/^i2)5 where M'l2 and 
are the dispersive and the absorptive part of the Bg — Bg mixing amplitudes, respectively. 
On the other hand, the parameters also contribute to the other processes, e.g. AMg, ATg 
and time- dependent CP asymmetries of the Bg meson decays. According to the global 
fit analysis in Ref. [1], the authors searched for the minimum for the processes which 
are sensitive to the Bg meson mixings as well as the like-sign dimuon charge asymmetry. 
Despite that the SM predicts a tiny angle giryM^^) ~ 10-(2-3), the result indicates 
an anomalously large CP phase of 0(0.1 — 1) in the Bg — Bs mixing [5]. The global fit 

^ The measurements of the Bg J /^^(j) decay have been recently updated by the Tevatron experi- 
ments [131 [H] • Although the results are turned to be consistent with the SM prediction, the experimental 
uncertainty is still large. 
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result also prefers a larger absorptive part of the Bs — Bg mixing than the SM prediction. 
However, both the experimental and SM values are largely uncertain [5l [16] . Thus, we 
focus on new physics (NP) contributions to the Bg — B^ mixing phase. 

In order to investigate the NP contribution to the Bg — Bg mixing, we parametrize 
dispersive part of the Bg — Bg mixing amplitude as 

M^, = (M^^)^^ + iMl,r^ = (M^2)^^(l + /^.e^-), (1) 

where (Mf2)^^ and (M^g)^^ the SM and NP contributions, respectively. This con- 
tributes to the mass difference of the neutral Bg mesons, AMg, as well as the semileptonic 
asymmetry, a^^. Since the measurements of AMg at the Tevatron are consistent with the 
SM prediction [I71[T8], (Mf2)'^^ is restricted as |1 + hge^^"""] ~ 1. On the other hand, ag is 
favored to be large by the measurement of the like-sign dimuon charge asymmetry. Thus, 
hg must be correlated with ag for ~ 1. In the absence of the NP contributions to the 
Bd — Bd mixing, the minimum is found around [5] 

a,) ^ (1.8,100°). (2) 

We will refer this value in the following discussion. Note that the following result does 
not change qualitatively as long as the experimental data favor a large mixing angle even 
if they shift in future. 

3 Gluino Dominant Case 

When superparticles have a mass of 0(100 — 1000)GeV, gluino-squark diagrams usually 
dominate SUSY contributions to the colored FCNC and CP violating processes. In the 
super-CKM basis [1], the down- and up-type squark mass matrices are described as 

(M|),, = diag(m;-),, + m%5%,, [Ml),, = dmg{ml)u + m2(5-),,, (3) 

where the first term in the right-hand side denotes the squark mass, and the second one 
represents the flavor-changing components with an average squark mass m?. Depending 
on the chirality structure of the squarks, the squark mixings are classified to the four 
types, (5x,i)ij, {.^'irdiji {^RL)ij i^RR)ijy with q = d,u and i,j the generation index. 
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Figure 1: Contours of hs = 2.0,1.5,1.0 and 0.5 from the outside to inside (solid line). 
The axes are the imaginary component of the b — s squark mixings. The inside of the 
dashed hne (circle) is allowed by Br(fe — )■ 57) at the 2a level. The region between the 
dotted-dashed (red) lines is the 2a band of Sri'x, and that inside the dotted (blue) lines 
is of S^K- We set msusY = = 500GeV, Re(5^i)23 = Re(4i?)23 = and tan/3 = 10. 

Then, the SUSY contributions to the CP-violating FCNCs are represented by the squark 
mixings and the masses of the squarks and gluino. It is noticed that the complex phase 
of the SUSY contributions is not rephasing invariant. We choose the phase convention 
such that the SM contributions are real, if not otherwise mentioned. Namely, the phase 
of represents the CP violation. 

The SUSY contribution to is dominated by the gluino box diagrams. Normalizing 
it by the SM contribution, /i^e^*""" is evaluated as [1] 

+«3 mn)l3 + (4l)23] - «4 [{S'MKM , (4) 

where the coefficients depend on the detail of the mass spectrum of the squarks and gluino, 
whose size is ai = 0(1), a2 = O(IOO), = 0(10) and 04 = 0(10) for rriq ~ nig ~ 500GeV 
(see e.g. for details). Let us mention that they are scaled as aj oc (5'^)23'^^ft) '^hile it is 
insensitive to tan/3 when the gluino diagram is dominant. On the other hand, it is noted 
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that the second hne of (jlj) consists of the squark mixings with chirahty flips. Since they 
are tightly constrained by Br(6 — 57) as {^j^n^Ril'i'i < O(10~^) for nig, nig ~ SOOGeV [1], 
their contributions are negligible. Thus, we assume {5lrrl)'2^ — simplicity. 

We show contours of hs in Fig. [H where we set rriq = nig = SOOGeV and Re((5^^)23 = 
^^i^RR)'23 — 0- is stressed that hg is enhanced by a combination of (5^^)23 and (5^jj)23 
as is found in (j4]). Thus, it becomes comparable to the SM contribution especially when 
the mixings are (5^^)23 ~ (<^^_r)23 ~ 0.1. For instance, we obtain hs = 1.8 for 1(5^^)231 = 
1(5^^)231 ^ 0.08 in the case of niq = mg = SOOGeV. 

As was mentioned in the previous section, the like-sign dimuon charge asymmetry 
favors a large Bg — Bg mixing phase. When the 02 term dominates the SUSY contributions 
in (jl]), the following angle gives the CP phase. 



where a, is related to 9+ as 6+ = 2(7, - 180°, e.g. 6+ ~ 20° for (x, = 100°. 

Before proceeding to the experimental constraints, let us comment on the SUSY contri- 
bution to r^2- Although the dimuon anomaly implies a larger than the SM value [51 [T6]. 
it is unlikely to expect this fulfilled in the SUSY models (see [T6]). 

The large squark mixings are restricted by the other observables. In order to obtain 
the minimum ([2]), the processes are considered which are sensitive to the new physics 
contributions to the Bg — Bg mixing. Moreover, the squark mixings (5"^) 23 can significantly 
contribute to the following observables: 

1. the branching ratio of the inclusive 6 — )■ 57 decay: Br(fe — )■ 57), 

2. the time-dependent CP asymmetry of the Bd decay into (pK and rj'K: S^Ks: Sn'Ks) 

3. the atomic electric dipole moments through the strange quark CEDM: dfg. 

Note that their experimental result agrees with the SM predictions. Thus, they can give 
additional constraint on (5^^)23 and (5^/j)23- 

Let us start from Br (6 — 57). The experimental result [19] is consistent with the SM 
prediction [20]. In fact, the difference is limited in 



e+ = arg(<5^^)23 + arg(4^)23, 



(5) 



0.3 X 10"^ < ABr(6 57) < 1.4 x 10 



-4 



(6) 
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at the 2a level. In terms of the effective Hamiltonian, the branching ratio is approximately 
given by jCy^p + |C*77p, where C^^ is the Wilson coefficient of the photonic magnetic 
operator, Oj^ = j^rrih^Sia'^'^ PRbi)F^^. It is known that the dominant SUSY contribution 
to Cj^ behaves as 

(C77)SUSY OC (5^^)23 tan /3 m J,, (7) 

where mgoft is a typical soft mass scale, and L ^ i? for Cy^ (see [21] for details). 

In Fig. m we display the region which is allowed by Br (6 — 37). We want to mention 
that the analysis bases on msoft = SOOGeV and tan /3 = 10. It is found that compared with 
the SUSY contributions to the — mixing (jl]), those to Br (6 — )■ 57) ([7]) is enhanced 
when tan /3 is large and/or soft masses are small. Thus, the b ^ s'-f constraint can exclude 
hs ^ 1 for large tan /3 and small soft masses. 

It is noticed that since the SM contributes mainly to the Wilson coefficients Ci (not Q), 
the decay rate is sensitive to extra contributions to the real component of Ci. There are 
briefly two contributions. First, the real component of (5^/^)23 can be finite with satisfying 
([2]), which induces KeC-j-y. Secondly, the chargino, neutralino and charged Higgs diagrams 
can also contribute. In particular, those with the CKM matrix as a source of the flavor 
mixing interfere with the SM contribution. In Fig. [H we take account of the chargino and 
charged Higgs diagrams with msoft = SOOGeV. They can change depending on the detailed 
mass spectrum of the superparticles, and we checked that there is a wide parameter region 
where ([2]) is fulfilled with the constraint from Br (6 — )■ 57) satisfied. 

Secondly, the time-dependent CP asymmetries of the B^ decays provide another con- 
straint especially on the imaginary component of the squark mixings. In particular, those 
of Bd —7- (f)K and Bd — t- rj'K have been measured well in the B-factories as [I9] 



0.20 < S^K < 0.88, 0.45 < Sr^'K < 0.73. (8) 

at the 2a level. If there is no NP contribution, it is considered that they become equal to 
the result of B^ J/ipK, S^k = 0.67 ± 0.2 [19]. In the SUSY models, they receive new 
contributions mainly through C^g and CsG? which are the Wilson coefficients of the chromo 



^ This constraint may be relaxed by a hierarchical squark mass spectrum between the first two and 
the third generations [TO] , 
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magnetic operators. Similarly to ([7]), the dominant SUSY contribution satisfies [2T] 

(C'sgJsusy oc (5^j23tan/3m-|, (9) 

and L ^ R for Csg- In the analysis, we use the generalized factorization method [221123] to 
evaluate the chromo-magnetic contribution to the decay amplitude, where the momentum 
which is transferred by the gluon of Osg is set ~ M^/2. 

It is stressed that 5*0/^ and Srj'K have different dependence on (5^^)23 and (5^^)23- Since 
the final state (pK has an odd parity, the decay amplitude depends on Ci + Ci, while noting 
that rj'K is parity even, the amplitude is proportional to Cj — C*j (see [21]). As a result, 
we obtain the allowed region in Fig. [H We notice that Sri'K indicates {S'Il)23 ~ ('^^,r)23, 
while S^K does not give a stringent bound because of the large experimental uncertainty. 

Compared with the SUSY contributions to the Bg — Bg mixing (jl]), those to (Csg) and 
(Csg) increase when tan/3 is large and/or superparticles are light (see (E])). Thus, large 
tan P and/or small rrisoft can close the window to realize /i^ ~ 1. It is also commented that 
they are less sensitive to the real component of Csg and Csg, because the real component 
of the decay amplitude is dominated by the SM contribution. 

As a result of the above constraints from b — > 57, S^k and Sr^'K as well as those 
sensitive to the anomalous Bg — Bg mixing, it seems that the SUSY contribution to the 
Bg — Bg mixing can be as large as the SM contribution with a large CP phase if the 
squark mixings are (5^^)23 ~ {^rr)23 ~ 0.1 for msoft = O(100)GeV and moderate tan/3. 
Moreover, focusing on the dependence on [Sf^ 101)23 and msoft of hg and the constraints, 
a larger superparticle mass alleviates the constraints. 

However, it is too early to conclude that the Bg — Bg mixing phase can be large with 
satisfying all the experimental constraints. The atomic and neutron EDMs are sensitive 
to the CP-violating phase of the b — s squark mixings through the CEDM of the strange 
quark [9]. Based on the hadronic calculation and the experimental results [21], the strange 
quark CEDM is bounded as (see [25] ) 

eld'gl < 1 X 10-^^ecm. (10) 

On the other hand, if both (5^j^)23 and (5^/j.)23 are finite, is dominated by the gluino- 



8 



10-23 
10-24 
10-25 
10-26 
10-27 

10-3 10-2 10-1 1 

arg(8RR) 

Figure 2: as a function of arg((5^j:j)23 with arg(5^^)23 = 0. The other parameters are 
rriq = rrig = /i = SOOGeV, 1(5^^)231 = |(^^i?)23| = 0.08 and tan/3 = 10. The red hne shows 
e\d1\ = 1.0 X 10-25ecm (see (^). 

squark diagrams, and we obtain [9]. 

e\dt\ (X |(5^j23(4«)23|sine-tan/3m-^ (11) 

where 0~ denotes the CP- violating phase, which is defined as 

e- = arg(5^j23 - arg(4^)23. (12) 

Ahhough G~ should also depend on the phase of the gluino mass and the fi parameter, 
they are tightly limited to be zero by the EDMs. Thus, we hereafter focus on the phase 
of (5ij23 and (5^^)23. 

The EDMs provide a very severe constraint on the phase G~. In Fig. [21 \dg\ is evaluated 
as a function of arg((5^^)23 with arg(5^j^)23 = 0, i.e. 6~ = — arg(5^^)23. Here, we choose 
the size of the squark mixings as 1(5^^)231 = \{^rr)23\ = 0.08, which realize hg = 1.8. 
Therefore, in order to fulfill ([2]) arg(5^^)23 must be correlated with arg(5^^)23 at the 
percent level. Even if we allow the CEDM bound up to 1 x 10~^^ecm due to a potentially 
large hadronic uncertainty, the phase is still required to tuned within the 10% level. Thus, 
we find that it is unlikely to expect a large CP asymmetry in the Bg meson mixing unless 
the phases satisfy ^ 6~ (see also [26]). 
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It is emphasized that the EDM constraint cannot be avoided even if we increase 
superparticle masses. This feature is contrasted to the constraints from b — 57, S(f,K 
and Srj'K- This is because the gluino-squark contribution to (fTTl) is proportional to 
('^'^)23"^soft' 'which is the same structure as hs (HI). This means that the EDMs provide 
more robust constraint on the Bs — Bg mixing phase than h — )■ 57, S^fix and Sr^iK-, and it 
is considered that the B^ — B^ mixing phase is difficult to be large. 

The strict EDM bound may be reconciled with the large Bs — Bs mixing phase by 
supposing cancellations among the SUSY contributions to the EDMs. In fact, there are 
additional contributions such as those from the neutralino and the chargino, which can 
interfere with the gluino contribution. Moreover, the SUSY contribution to the Bs — Bs 
mixing phase, G"*", in ([5]) is not equal to the CP- violating phase appearing in the CEDM, 
6~, in (fT2!) . Thus, the CEDM bound can be avoided if 6~ vanishes, i.e. arg((5^^)23 — 
arg(5^j:j)23, while the Bs — Bs mixing phase is left sizable as long as arg(5^^)23 and 
arg((5^^)23 are large [lO]. 

This (accidental) cancellation is, however, easily polluted by radiative corrections. 
One of the most ubiquitous corrections is renormalization group contributions to the soft 
parameters. Even if we suppose (5^^)23 = (^rj?)23 at a high energy scale, (5^^)23 receives 
a correction during the evolution down to the weak scale as 



AKJ.3-iK:K,^>"^. (13) 



where Mx and Msusy are the input and SUSY scales, respectively. For instance, if 
Mx is the GUT scale, A {^'11)2^, becomes O(10~^). It is emphasized that the phase 
of the correction is aligned to the CKM phase, that is real in this article. Thus, this 



'23 I 



generally shifts the total phase of (5^j^)23 by 0(0.1 — 1) when the mixing is 1(5^^ 
(0.01 — 0.1). This means that it is quite unlikely to expect the cancellation among the 
SUSY contributions to the CEDM unless it is protected by symmetries. |^ Consequently, 
it is very difficult to reconcile the large Bs — Bs mixing phase with the EDM bound when 
the gluino diagrams dominate the SUSY contributions. 



■^When the CP violation originates from the Yukawa couphngs apart from the CKM matrix, the 
hermicity of the Yukawa matrix ensures the suppression of the EDMs in relation to the strong CP 
problem |10) . 
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4 Decoupling Case 



When the gluino is relatively light, we showed in the previous section that the large CP 
asymmetry in the Bs — mixing conflicts with the EDM bound generically. In order 
to avoid the constraint, let us consider heavy superparticle scenarios. When the colored 
superparticles are sufficiently heavy, the gluino-squark contributions are suppressed. 

There are possibly sizable contributions from the heavy Higgs bosons as long as their 
masses are 0(100 — 1000)GeV. It is stressed that the Higgs contributions are enhanced 
when tan/3 is large. The Lagrangian of the Higgs bosons and the down-type quarks is 
written as [121 US] 



where Yi is the Yukawa coupling of the down-type quark, i, and Y-j is a non-holomorphic 
Yukawa coupling, which is generally complex and flavor non-diagonal. Diagonalizing the 
down-type quark mass matrix, we obtain the flavor- changing (and CP-violating) neutral 
Higgs couplings as [121 CHI [H] 



It is noticed that the Higgs bosons couple to the down-type quarks with flavor violations 
due to the non-holomorphic interaction. 

Within the framework of SUSY, the non-holomorphic interactions are not involved in 
the Lagrangian at the tree level. Rather, they are generated by radiative corrections to the 
Higgs coupling with the matter fermions, whose diagrams include the superparticles [121 
[T3] . The important features of them are that the non-holomorphic Yukawa couplings 
are proportional to tan /3 and are not suppressed when the superparticles decouple. As a 
result, the second term in ( !T5|l is proportional to tan^ /3 even in the decoupling limit of the 
superparticles. It is also noted that the heavy Higgs bosons mainly consist of when 
their masses are relatively large compared to the electroweak scale. Thus, the diagrams 
intermediating the heavy Higgs bosons are sizably enhanced when tan f3 is large, even 
though the couplings emerges at the one-loop order. 

The neutral heavy Higgs diagram can enhance the Bg — Bg mixing. The double Higgs 
penguin diagrams dominantly contribute to which include a couple of the heavy Higgs 



C = Y,Q,H,D'r + Y^Q^KD'^ + h.c. 



(14) 



£ = YljQiHlD'-j - Yl^ tan /3 QiD^^Hd + h.c. 



(15) 
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Figure 3: Contours of hs = 2.0, 1.5, 1.0 and 0.5 from the outside to inside. The constraint 
from Br (5 — fi~^fi~) is shown by the dashed (red) hne. We take rriH = 500GeV and 
Re(5^^) = Re(4^) = 0. In the left (right) panel, we set tan/3 = 30(40). 




0.05 0.1 0.15 0.2 



Im(5iJ=Im(4^) 

Figure 4: Bt{Bs — )■ /i'''/i~) for tan/3 = 30 and 40, where the squark mixings are varied 
with Im((i^^)23 = Ini('^R_R)23 fixed. The parameters are Re((5^^)23 = Re((5^^)23 = 
and rriH = 500GeV. The dash-dotted (red) line shows the current experimental bound, 
Bt{Bs ;U+/i-) = 4.7 X 10"^ [27], and the dotted-line (blue) is the SM prediction. 
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couplings. Then, Mfj is roughly enhanced by tan^ /3. The contribution becomes [28] l 



a[V:MR)23 - 4(5^l)23(5^«)23 , (16) 



where we explicitly show the flavor mixings which come from the non-holomorphic cou- 
plings. In particular, the flrst term in the right-hand side consists of the chargino and 
gluino contributions to the non-holomorphic couplings, and the second term originates 
from the gluinos. The coefficients a[ include the loop functions, and they are estimated 
to be a real and positive number. To be explicit, they depend on parameters as 



tan^ f3 
[I + etan/3)4m2^' 



0^ T^^Tirml— 2-' (1^) 



where the heavy Higgs mass, mn, arises from the contributions from the CP-even and 
CP-odd heavy Higgs bosons. Quantitatively, it is estimated as a[ ~ a'g. In contrast 
to the gluino dominant case, we notice that hg can be sizable solely by (5^/^)23, while 
the contribution which depends only on (5^^)23 is suppressed by a chiral factor ms/rrih. 
Therefore, a large hs favores large (5r/j)23- 

The constant e comes from the non-holomorphic corrections to the Yukawa couplings 
of the down-type quarks. It satisfies Y-'- ~ y^e, and e is not enhanced by tan /3 (see [28] 
for details). When superparticles are almost degenerate, it is estimated as |e| ~ 0.01, and 
its sign is determined by the sign of fi as sgn(e) = sgn(/i). As we will explain later, since 
negative /i leads to larger d'^ and Br (6 — )■ 57), positive sign is preferred. In this article we 
take /i > 0, if not otherwise mentioned. 

We shows contours of in Fig. [3l where we take Re(5^^)23 = Re{6'^j^)23 = and 
rriH = SOOGeV. We also choose two cases of tan /3 as tan /3 = 30 and 40 in the left and 
right panels, respectively. It is found that the heavy Higgs contributions to the Bs — Bg 
mixing become comparable to the SM contribution for {Sfi jip/) 23 ~ 0.1. For example, we 
obtain hs = 1.8 when the mixings are Im((5^^)23 — 0.07(0.03) for tan/3 = 30(40). 

The fiavor-changing neutral Higgs couplings also enhance the branching ratio of the 
Bg — fi^fJ'" decay, significantly. The dominant contribution to the decay amplitude 
includes the non-holomorphic interactions and the Yukawa coupling of the muon. Thus, 



^In the numerical analysis, we also include the charged Higgs contributions [2H], though they are 
subdominant compared to the neutral Higgs contributions. 
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Br(i?s — !■ fi~^fi ) increases at tan^/3. In terms of the effective Hamiltonian, tlie branching 
ratio is represented by the Wilson coefficients as |28] 

(18) 

where Asm denotes the SM contribution. The Wilson coefficients Cs and Cs are induced 
by the heavy Higgs bosons and given as 

In the bracket, we explicitly show the flavor mixings and a relative size of the contributions 
by the coefficients, a and b, which are estimated as a ~ 0.1 and 6 ~ 1, respectively. 
Numerically, Cs and Cs exceed Asm when the squark mixings are {Sf^ rb)23 ~ 0.1. Then, 
Bi{Bs is roughly proportional to 1(5^^)23^ + |(^fl'R)23p- 

The Bs — )■ fJ'~^fJ'~ decay has not been measured yet, and the branching ratio is con- 
strained by the Tevatron experiments as 

Br{Bs < 4.7 X 10"^ (20) 

at the 90% C.L. [27]. In Fig. [3l we show the constraint of Bt{Bs — > /i~^/i~). Namely, the 
region inside the ellipse is consistent with the experiment. We notice that the contribution 
proportional to (5^^)23 becomes smaller than those to (5^^)23? because the sign of a is 
opposite to that of b in Cs- Therefore, the allowed region for (5^/,)23 is rather large in 
Fig. O Anyway, we obtain a large Bg — Bg mixing without conflicting with the constraint 
of Bt{Bs —7- fJ''^fJ'~) in a wide parameter region. 

In the LHCb experiment, it is expected that the sensitivity of Bt{Bs — t- fi^fi~) reaches 
the SM prediction [29], which is Br(5, j^^j^^) = (3.35 ± 0.32) x 10^^ [30]. In Fig. H 
Bt{Bs —7- /i"*'/i~) is shown as a function of the size of the mixing with Im(5^j^)23 = 
Im(5^j:j)23 fixed. When {^'Il^rr)^?, is large enough to satisfy ([2]), the branching ratio is 
found to be close to the current upper bound for tan /3 = 40. It is noted that Br(i?s — )■ 
/i"*"/i~) depends on the heavy Higgs boson mass. We set the mass to be SOOGeV in 
Fig. m Compared with the dependence of hs on mn and ((5'^)23, Bi{Bs fJ'~^fJ'~) relatively 
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Figure 5: as a function of arg((5^^)23. We take arg(5^^)23 = and mH± = SOOGeV. The 
solid (dashed) hue corresponds to 1(5^^)231 = |(^^r)23| = 0.07(0.03) for tan/3 = 40(30). 
The dash-dotted (red) hne shows e|(i^| = 1.0 x 10~^^ecm. 

decreases as rriH becomes larger with hs fixed. Nonetheless, Bt{Bs /i'^/i") tends to be 
sizable to obtain a large 5^ — i^o^mixing phase, and such a large branching ratio is expected 
to be measured in the LHCb. □ 

Let us turn to the experimental constraints from the EDMs, Br (6 — >■ 57), Srj'Ks ^ind 
S^Ks- the gluino dominant case, the EDMs provide especially tight constraints on 
the CP violations. In contrast, when the superparticles are heavy, we will see that these 
constraints are greatly ameliorated. 

In the decoupling limit, these processes are induced by the charged Higgs diagrams. 
Similarly to the heavy neutral Higgs bosons, the flavor- changing (CP- violating) charged 
Higgs contributions are generated by the non-holomorphic interactions. Since the mass 
of the charged Higgs boson, mH±, is likely to degenerate with that of the heavy Higgs 
bosons up to SU(2) breaking corrections, we hereafter assume mf^± = rriH for simplicity. 

The charged Higgs contribution to the CEDM of the strange quark is given at the 
BLO as [31]. 

where the CKM matrix gives arise at the charged Higgs vertex, and the squark mix- 

^Note that since the contributions from the real part of {6'^)23 interferes with SM and minimal flavor- 
violating contributions, Br(i?s /^+/i^) can be affected by the real component. 
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ing originates in the the non-holomorphic Yukawa couphng which is dominated by the 
gluino diagrams. In addition, there is a contribution which depends on [{S'[i,)23{^'rr)23]- 
This term is obtained by substituting the vertex of the CKM matrix into another non- 
holomorphic couphng. Since the non-holomorphic interaction is smaller than the top 
Yukawa coupling, its contribution to is smaller by one order of magnitude than (|2T]1 . 
Thus, the EDM is sensitive to the phase of (5^/j)23- Although the neutral Higgs bosons 
can also contribute to d^, the effects are negligible because the correction from the CP- 
even heavy Higgs boson cancels with that from the CP-odd heavy Higgs boson when the 
mass splitting between rriA and tuh is small [3T] . 

Fig. Elshows the dependence of d'^ on the phase of ((5^i?)23 in the cases of tan (3 = 30 and 
40. In the analysis, the mixings are chosen as (^^75)23 = 0.07 and 0.03, respectively, which 
are required to realize hg = 1.8. It is emphasized that d'^ is proportional to (5'')23 tan /3. 
This is contrasted to the SUSY contributions to the B — B mixing and Br(i?s — t- 
which briefly depend on (5'')23 tan^ (5 and tan^ /3, respectively. Moreover, d1 is scaled 
by which is the same as hs- Thus, when hg is given, the EDM bound is satisfied 

especially when tan /3 and 'mu± are larger. Actually, for tan /3 = 40 and {5jifi)2z = 0.03, 
e|(i^| is smaller than 1.0 x 10~^^ecm even with arg(5^j:j)23 = 0(1). Then, Br(i?5 — 
is likely to be large in the region. 

The charged Higgs contribution to c?^ is enhanced when e is negative due to the factor 
(l-|-etan/3)~^ in f l^T]) . As is the case of /i^, it comes from the non-holomorphic corrections 
to the Yukawa couplings of the down-type quarks. Note that the sign of e is equal to that 
of /i. Hence, when /i is negative, c?* becomes larger by 0(10) than the result of positive 
/i. As a result, can exceed the experimental bound for axg{5'^pf)22, = 0(1). Thus, /i > 
is preferred by the EDMs. 

The charged Higgs boson also contributes to Br (6 — 57), Sr^'Xg and S^Ks through Cj^ 
and CsG- At the BLO, Oy^ and Csg become 



l-6tan/3 1 ,.^(.r.-u^_^aii^_(j:d^ 1 
^^-^^ l + etan/3<: + ^('° ^ ^ 1 + etan/3 ^^^^^^^ 

tan/3 ^ 1 
l^Tt^^^^^^^'^^i 



where the first term in the right-hand side of C^^ sg comes from the charged Higgs diagram 
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with the CKM matrix. On the other hand, the flavor mixing in the second term and that 
of Ci^f^sG originate from the non-holomorphic interaction which is dominated by the gluino 
diagrams. 

Ahhough we omit the numerical coefficients, it is checked that (l22l) is quantitatively 
smaller than the SM contribution for fi > 0. Actually, the first term of Cj^^sg becomes 
O(10~^) for mH± = SOOGeV and tan/3 = 30 — 40, which is just a small correction to the 
SM value, {C7.y^sG)sM ~ 0.1. On the other hand, the other terms are linearly proportional 
to tan /3, in contrast to hg and Bt{Bs — )■ fi^fi~). Numerically, we obtain Cj^^sG ~ 10~^ 
for tan/3 = 0(10) and (5'^)23 ~ 0.1 from these terms. As a result, Br(6 — )■ 57), S^Ks 
Sn'Ks remain consistent with the experimental values. 

When /i is negative, Cj^^sc becomes larger that the case of /i > 0. Actually, they can 
reach O{10~^) for large tan /3 mainly from the first term in ( 122|) . Since this is comparable to 
the SM values, the charged Higgs contributions to Br (6 — > 57) can exceed the experimental 
bound. On the other hand, noting that the first term in f l22|) is real, i.e. aligned to the 
SM contribution, this contribution to Stj^Ks ^i-iid S^k^ is irrelevant. It is also checked that 
the other contributions in f l22p are less significant for Srj'Ks S^Ks^ because they are 
smaller than the first term. 

5 Conclusion 

Motivated by the like-sign dimuon charge asymmetry, we have investigated the CP asym- 
metry in the Bs — Bg mixing. Although a large mixing phase is consistent with the 
experimental constraints from Br(6 — 57) and the mixing-induced CP asymmetries of 
the Bd decays into (pK and t]'K, we have shown that when the gluino diagrams domi- 
nate the SUSY contributions, the EDM bound almost excludes this possibility unless the 
cancellation works among the SUSY contributions to the EDM. 

In order to alleviate the tight EDM constraint, we have analyzed the heavy superpar- 
ticle scenarios. In the decoupling limit of the superparticles except for the Higgs bosons, 
the FCNCs and CP violations are indued by intermediating the neutral/charged heavy 
Higgs bosons. We have found that the non-holomorphic Yukawa couplings can enhance 
the Bg — Bg mixing with satisfying the experimental constraints from the EDMs as well 
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as Br (6 — )■ 57), S^k and Sr^iK- In the region where the Bg — Bg mixing phase is large, 
Br(i?s — tends to be sizable, and we expect to observe the decay at the LHCb. 
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